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Abstract 

This is the extended version of the preprint based on the lectures given in 
Cargese Summer School and Chernogolovka Summer School in 93. The incompress- 
ible fluid dynamics is reformulated as dynamics of closed loops C in coordinate 
space. We derive explicit functional equation for the pdf of the circulation PciX) 
which allows the scaling solutions in inertial range of spatial scales. The pdf decays 
as exponential of some power of /A^ where A is the minimal area inside the loop. 
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1 Introduction 



Incompressible fluid dynamics underlies the vast majority of natural phenomena. It is 
described by famous Navier-Stokes equation 

= vdlvc - vpdpVa - daP; daVa = (1) 

which is nonlinear, and therefore hard to solve. This nonlinearity makes life more inter- 
esting, though, as it leads to turbulence. Solving this equation with appropriate initial 
and boundary conditions we expect to obtain the chaotic behavior of velocity field. 

The simplest boundary conditions correspond to infinite space with vanishing velocity 
at infinity. We are looking for the translation invariant probability distribution for velocity 
field, with infinite range of the wavelengths. In order to compensate for the energy 
dissipation, we add the usual random force to the Navier-Stokes equations, with the short 
wavelength support, corresponding to large scale energy pumping. 

One may attempt to describe this probability distribution by the Hopf generating func- 
tional (the angular bracket denote time averaging, or ensemble averaging over realizations 
of the random forces) 

Z[J] = ^exp Q (frJ^{r)v^{r)Yj (2) 
which is known to satisfy linear functional differential equation 

Z (3) 

similar to the Schrodinger equation for Quantum Field Theory, and equally hard to solve. 
Nobody managed to go beyond the Taylor expansion in source J , which corresponds to 
the obvious chain of equations for the equal time correlation functions of velocity field in 
various points in space. The same equations could be obtained directly from Navier-Stokes 
equations, so the Hopf equation looks useless. 

In this workQ we argue, that one could significantly simplify the Hopf functional with- 
out loosing information about correlation functions. This simplified functional depends 
upon the set of 3 periodic functions of one variable 

C : r„ = Ca(^); < ^ < 27r (4) 

which set describes the closed loop in coordinate space. The correlation functions reduce 
to certain functional derivatives of our loop functional with respect to C{9) at vanishing 
loop C ^ 0. 

The properties of the loop functional at large loop C also have physical significance. 
Like the Wilson loops in Gauge Theory, they describe the statistics of large scale structures 
of vorticity field, which is analogous to the gauge field strength. 

^see also where this approach was initiated and ||] where its relation with the generalized Hamil- 
tonian dynamics and the Gibbs-Boltzmann statistics was established 
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In Appendix A we recover the expansion in inverse powers of viscosity by direct iter- 
ations of the loop equation. 

In Appendix B we study the matrix formulation of the Navier-Stokes equation, which 
may serve as a basis of the random matrix description of turbulence. 

In Appendix C we study the reduced dynamics, corresponding to the functional Fourier 
transform of the loop functional. We argue, that instead of 3D Navier-Stokes equations 
one can use the ID equations for the Fourier loop Pa{0, t). 

In Appendix D we discuss the relation between the initial data for velocity field and 
the P field, and we find particular realisation for these initial data in terms of the gaussian 
random variables. 

In Appendix E we introduce the generating functional for the scalar products Pa{O)Pa{0'). 
The advantage of this functional over the original "^[C] functional is the smoother con- 
tinuum limit. 

In Appendix F we discuss the possible numerical implementations of the reduced loop 
dynamics. 

In Appendix G we show uniqueness of the tensor area law within certain class of 
functionals. 

In Appendix H we present the modern view at the old problem of the minimal surface. 

In Appendix I we show that the triple Kolmogorov correlation function corresponds 
to a vanishing correlation of vorticity with two velocity fields. 

2 The Loop Calculus 

We suggest to use in turbulence the following version of the Hopf functional 

vi>[c] = (exp(^ j dc^{e)v^{c{e))yj (5) 

which we call the loop functional or the loop field. It is implied that all angular variable 6 
run from to 27rand that all the functions of this variable are 27r periodic. The viscosity 
V was inserted in denominator in exponential, as the only parameter of proper dimension. 
As we shall see below, it plays the role, similar to the Planck's constant in Quantum 
mechanics, the turbulencecorresponding to the WKB limit — 0. Q 

As for the imaginary unit i , there are two reasons to insert it in the exponential. 
First, it makes the motion compact: the phase factor goes around the unit circle, when 
the velocity field fluctuates. So, at large times one may expect the ergodicity, with well 

^This parametrization of the loop is a matter of convention, as the loop functional is parametric 
invariant. 

•^One could also insert any numerical parameter in exponential, but this factor could be eliminated by 
space- and/or time rescaling. 
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defined average functional bounded by 1 by absolute value. Second, with this factor of 
t , the irreversibility of the problem is manifest. The time reversal corresponds to the 
complex conjugation of so that imaginary part of the asymptotic value of \1/ at t oo 
measures the effects of dissipation. 

The loop orientation reversal C{9) C{2n — 9) also leads to the complex conjugation, 
so it is equivalent to the time reversal. This symmetry implies, that any correlator of 
odd/even number of velocities should be integrated odd/even number of times over the 
loop, and it must enter with an imaginary /real factor. Later, we shall use this property 
in the area law. 

We shall often use the field theory notations for the loop integrals, 

^[q = ^exp(^^£dr,«,)^ (6) 
This loop integral can be reduced to the surface integral of vorticity field 

i^f,u = df^Vu - dyV^, (7) 

by the Stokes theorem 

TcH = <j>^ draVa = da^i^uj^y] dS = C (8) 

This is the well-known velocity circulation, which measures the net strength of the 
vortex lines, passing through the loop C. Would we fix initial loop C and let it move 
with the flow, the loop field would be conserved by the Euler equation, so that only the 
viscosity effects would be responsible for its time evolution. However, this is not what 
we are trying to do. We take the Euler rather than Lagrange dynamics, so that the loop 
is fixed in space, and hence \l/ is time dependent already in the Euler equations. The 
difference between Euler and Navier-Stokes equations is the time irreversibility, which 
leads to complex average \l/ in Navier-Stokes dynamics. 

It is implied that this field \l/ [C] is invariant under translations of the loop C{9) — > 
C{9) + const. The asymptotic behavior at large time with proper random forcing reaches 
certain fixed point, governed by the translation- and scale invariant equations, which we 
derive in this paper. 

The general Hopf functional (0) reduces for the loop field for the following imaginary 
singular source 

Ja{r) = -i dr'j' (r' - r) (9) 
v Jc 

The \l/ functional involves connected correlation functions of the powers of circulation 
at equal times. 

^[C]=exp(|:-^((r"^H))) (10) 
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This expansion goes in powers of effective Reynolds number, so it diverges in turbulent 
region. There, the opposite WKB approximation will be used. 

Let us come back to the general case of the arbitrary Reynolds number. What could 
be the use of such restricted Hopf functional? At first glance it seems that we lost most 
of information, described by the Hopf functional, as the general Hopf source J depends 
upon 3 variables x, y, z whereas the loop C depends of only one parameter 6. Still, this 
information can be recovered by taking the loops of the singular shape, such as two 
infinitesimal loops Ri,R2, connected by a couple of wires 

G>^ >© 

Fig. 1 



The loop field in this case reduces to 



* [C] ^ (exp ( _Ej>,.(rO + -S«>,.(r2) ) > (11) 



where 



Sj. = / dr^r, (12) 

is the tensor area inside the loop R. Taking functional derivatives with respect to the 
shape of Ri and R2 prior to shrinking them to points, we can bring down the product of 
vorticities at ri and r2. Namely, the variations yield 

5S^^ = i {drjr^ + r^dSr^) = i {dr„5r^ - drf,5r„) (13) 

•J R -J Fi 

where integration by parts was used in the second term. 

One may introduce the area derivative j^^-jp^i which brings down the vorticity at the 
given point r at the loop. 

The careful definition of these area derivatives are or paramount importance to us. 
The corresponding loop calculus was developed inP] in the context of the gauge theory. 
Here we rephrase and further refine the definitions and relations established in that work. 

The basic element of the loop calculus is what we suggest to call the spike derivative, 
namely the operator which adds the infinitesimal A shaped spike to the loop 

_ (15) 
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The finite spike operator 

A(r,^,e) = exp(raD„(^,e)) (16) 

adds the spike of the height r. This is the straight hne from C{9) to C{9 + €)+r, followed 
by another straight hne from C(^ + e)+r to C(^ + 2e), 




Fig. 2 

Note, that the loop remains closed, and the slopes remain finite, only the second 
derivatives diverge. The continuity and closure of the loop eliminates the potential part 
of velocity; as we shall see below, this is necessary to obtain the loop equation. 

In the hmit e — > these spikes are invisible, at least for the smooth vorticity field, as 
one can see from the Stokes theorem (the area inside the spike goes to zero as e) . However, 
taking certain derivatives prior to the limit e — > we can obtain the finite contribution. 

Let us consider the operator 

n(r,r',^,e)=A(r,^,^e)A(r',^,e) (17) 

By construction it inserts the smaller spike on top of a bigger one, in such a way, that a 
polygon appears 




Fig. 3 

Taking the derivatives with respect to the vertices of this polygon r, r' , setting r ~ 
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r' = and antisymmetrising, we find the tensor operator 

^o.p{e, e) = -I uD^ (e, ie) Dfs (9, e) - {a ^ /?} (18) 
which brings down the vorticity, when applied to the loop field 

n^p{e, e)^ [C] ^ ujo^p {cm ^ [C] (19) 

The quick way to check these formulas is to use formal functional derivatives 

Taking one more functional derivative derivative we find the term with vorticity times 
first derivative of the 6 function, coming from the variation of C'{9) 



SC^{9)6Ct,{9') ' '6a^p{C{9)) ' ' Sa^^ {C{9)) 6ap^{C{9')) 

This term is the only one, which survives the limit e ^ in our relation (p!9|). 

So, the area derivative can be defined from the antisymmetric tensor part of the second 
functional derivative as the coefficient in front of 6'{9 — 9') . Still, it has all the properties 
of the first functional derivative, as it can also be defined from the above first variation. 
The advantage of dealing with spikes is the control over the limit e — > , which might be 
quite singular in applications. 

So far we managed to insert the vorticity at the loop C by variations of the loop field. 
Later we shall need the vorticity off the loop, in arbitrary point in space. This can be 
achieved by the following combination of the spike operators 

A{r,9,e)U{n,r,,9 + e,6); 5«e (22) 

This operator inserts the 11 shaped little loop at the top of the bigger spike, in other 
words, this little loop is translated by a distance r by the big spike. 

Taking derivatives, we find the operator of finite translation of the vorticity 

A{r,9,e)Q^p{9 + e,d) (23) 
and the corresponding infinitesimal translation operator 

D^{9,e)n^^{9 + e,6) (24) 

which inserts d^uOap {C{9)) when apphed to the loop field. 

Coming back to the correlation function, we are going now to construct the operator, 
which would insert two vorticities separated by a distance. Let us note that the global A 
spike 

A (,-. 0, = e.p L ^) ^) (25) 



when applied to a shrunk loop C((/>) = does nothing but the backtracking from to r 

< > r 

Fig. 4 

This means that the operator 

f^a/3(0,5)l^Ap(vr,5)A(r,0,7r) (26) 
when applied to the loop field for a shrunk loop yields the vorticity correlation function 

^ap{0, S)n^p{7r, 6)A (r, 0, tt) ^[0] = MOpUr)) (27) 

The higher correlation functions of vorticities could be constructed in a similar fashion, 
using the spike operators. As for the velocity, one should solve the Poisson equation 



with the proper boundary conditions , say, v — at infinity. Formally, 

This suggests the following formal definition of the velocity operator 

1 



(28) 



(29) 



Another version of this formula is the following integral 

V^{9, e,S)=J (^V^^A (p, 9, e) n^piO + e, S) 



(30) 
(31) 

(32) 



where the A operator shifts the Q by a distance p off the original loop at the point 
r^C{e + e) 




Fig. 5 
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3 Loop Equation 



Let us now derive exact equation for the loop functional. Taking the time derivative of 
the original definition, and using the Navier-Stokes equation we get in front of exponential 

dra^ (vdlv^ - vpdpVa - dap) (33) 

The term with the pressure gradient yields zero after integration over the closed loop, and 
the velocity gradients in the first two terms could be expressed in terms of vorticity up 
to irrelevant gradient terms, so that we find 

<j)^ drj- {udpupa - vpUfSa) (34) 

Replacing the vorticity and velocity by the operators discussed in the previous Section 
we find the following loop equation (in explicit notations) 

-i^[C]= (35) 

The more compact form of this equation, using the notations of , reads 

I z/^[C] = T^c^ (36) 

Jc \ oapair) J 4:7r\r - r'\-^ daf3a{r)daf3.^{r') J 



Now we observe that viscosity u appears in front of time and spatial derivatives, like the 
Planck constant h in Quantum mechanics. Our loop hamiltonian Tic is not hermitean, 
due to dissipation. It contains the second loop derivatives, so it represents a (nonlocal!) 
kinetic term in loop space. 

So far, we considered so called decaying turbulence, without external energy source. 
The energy 

E = J dVi vl (37) 

would eventually all dissipate, so that the fluid would stop. In this case the loop wave 
function would asymptotically approach 1 

^C] 1 (38) 

In order to reach the steady state, we add to the right side of the Navier-Stokes 
equation the usual gaussian random forces /«(?", t) with the space dependent correlation 
function 

{fair, t)ff,{r', t')) = 5^p5{t - t')F{r - r') (39) 
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concentrated at at small wavelengths, i.e. slowly varying with r — r' . 
Using the identity 



(/„(r,t)$[t;(.)]) = j d?r'F{T-T')^-^^ (40) 

which is valid for arbitrary functional $ we find the following imaginary potential term 
in the loop hamiltonian 

SHc = I U[C] = - (f dr^i dr'^F{r - r') (41) 

Note, that orientation reversal together with complex conjugation changes the sign of 
the loop hamiltonian, as it should. The potential part involves two loop integrations times 
imaginary constant. The first term in the kinetic part has one loop integration, one loop 
derivative times imaginary constant. The second kinetic term has one loop integration, 
two loop derivatives and real constant. The left side of the loop equation has no loop 
integrations, no loop derivatives, but has a factor of t . 

The relation between the potential and kinetic parts of the loop hamiltonian depends 
of viscosity, or, better to say, it depends upon the Reynolds number, which is the ratio 
of the typical circulation to viscosity. In the viscous limit, when the Reynolds number 
is small, the loop wave function is close to 1. The perturbation expansion in ^ goes in 
powers of the potential, in the same way, as in Quantum mechanics. The second (nonlocal) 
term in kinetic part of the hamiltonian also serves as a small perturbation (it corresponds 
to nonlinear term in the Navier-Stokes equation). The first term of this perturbation 
expansion is just 



d^k F{k) 



<j> dVaC' 



kr 



(42) 



(27r)3 2z/3 P 

with F{k) being the Fourier transform of F{r). This term is real, as it corresponds to the 
two- velocity correlation. The next term comes from the triple correlation of velocity, and 
this term is purely imaginary, so that the dissipation shows up. 

This expansion can be derived by direct iterations in the loop space as in [0], inverting 
the operator in the local part of the kinetic term in the hamiltonian. This expansion 
is discussed in Appendix A. The results agree with the straightforward iterations of the 
Navier-Stokes equations in powers of the random force, starting from zero velocity. 

So, we have the familiar situation, like in QCD, where the perturbation theory breaks 
because of the infrared divergencies. For arbitrarily small force, in a large system, the 
region of small k would yield large contribution to the terms of the perturbation expansion. 
Therefore, one should take the opposite WKB limit — >^ 0. 

In this limit, the wave function should behave as the usual WKB wave function, i.e. 
as an exponential 

vl/[q^exp(i^) (43) 
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The effective loop Action S[C] satisfies the loop space Hamilton- Jacobi equation 

m ^ -^Un + £,.„/,3,._^_^_i|_ (44) 

The imaginary part of S[C] comes from imaginary potential U[C], which distinguishes 
our theory from the reversible Quantum mechanics. The sign of 535" must be positive 
definite, since < 1. As for the real part of S[C], it changes the sign under the loop 
orientation reversal C{9) C{2tx — 6). 

At finite viscosity there would be an additional term 

- V i dr^ds^^ -ivldrJ (45) 

JC " ^o-fjair) Jc J 47r|r — r'P 5(T/3Q,(r)5(T^^(r') 

on the right of (0). As for the term 

- I dr^i {dpS[C]) (46) 

which formally arises in the loop equation, this term vanishes, since dj3S[C] = 0. This 
operator inserts backtracking at some point at the loop without first applying the loop 
derivative at this point. As it was discussed in the previous Section, such backtracking 
does not change the loop functional. This issue was discussed at length in where the 
Leibnitz rule for the operator daj^ was established 

a « . (47) 

In other words, this operator acts as a first order derivative on the loop functional with 
finite area derivative (so called Stokes type functional). Then, the above term does not 
appear. 

The Action functional S[C] describes the distribution of the large scale vorticity struc- 
tures, and hence it should not depend of viscosity. In terms of the above connected 
correlation functions of the circulation this corresponds to the limit, when the effective 
Reynolds number ^^^^ goes to infinity, but the sum of the divergent series tends to the 
finite limit. According to the standard picture of turbulence, the large scale vorticity 
structures depend upon the energy pumping, rather than the energy dissipation. 

It is understood that both time t and the loop sizeQ \C\ should be greater then the 
viscous scales 

t > to = iy^£'^] \C\ > ro = z/t^-i (48) 
where £ is the energy dissipation rate. 

It is defined from the energy balance equation 

= dt(^^vlj = u {v^dW) + {f^v^) (49) 



*As a measure of the loop size one may take the square root of the minimal area inside the loop. 
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which can be transformed to ^ 

-v{ui;) = ^F{Q) (50) 

The left side represents the energy, dissipated at small scale due to viscosity, and the 
right side - the energy pumped in from the large scales due to the random forces. Their 
common value is £. 

We see, that constant F{r — r'), i.e., F(k) oc S{k) is sufficient to provide the necessary 
energy pumping. However, such forcing does not produce vorticity, which we readily see 
in our equation. The contribution from this constant part to the potential in our loop 
equation drops out (this is a closed loop integral of total derivative). This is important, 
because this term would have the wrong order of magnitude in the turbulent limit - it 
would grow as the Reynolds number. 

Dropping this term, we arrive at remarkably simple and universal functional equation 

The stationary solution of this equation describes the steady distribution of the circulation 
in the strong turbulence. Note, that the stationary solutions come in pairs ±5*. The sign 
should be chosen so, that SS* > 0, to provide the inequality l^'l < 1. 



4 Scaling law 

The 'Hamilton- Jacobi' equation without the potential term 
scaling solutions 

S[C] = t^^-^c 



c 



allows the family of the 
(52) 



with arbitrary index n. The scaling function satisfies the equation 



Jc dapa[r) Jc 



dr^ / d r 



6(f)[C] <50[C] 



(53) 



47r|r — r'P 5cr/3Q,(r) 5af3^{r' 

The left side here was computed, using the chain rule differentiation of functional. 

Asymptotically, at large time, we expect the fixed point, which is the homogeneous 
functional 

SooiC] = \c\'-^f 

zeroing the right side of our 'kinetic' functional equation 

6S^[C] 6S^[C] 



C_ 



(54) 







dr^ 



c 



d\' ""^ 



47r|r - r'P 5(ypa{r) 5ap^{r') 



The Kolmogorov scaling would correspond to 

3 

K = - 

2 



(55) 



(56) 
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in which case one can express the S functional in terms of S 

s[c] = se' ^ 



(57) 



One can easily rephrase the Kolmogorov arguments in the loop space. The relation 
between the energy dissipation rate and the velocity correlator reads 

£ = {va{ro)vp{0)dpVam (58) 

where the point splitting at the viscous scale tq is introduced. Such splitting is necessary 
to avoid the viscosity effects; without the splitting the average would formally reduce to 
the total derivative and vanish. 

Instead of the point splitting one may introduce the finite loop of the viscous scale 
\C\ ~ To, and compute this correlator in presence of such loop. This reduces to the WKB 
estimates 

u^a,ir) - v^{r) = f d\' /^~\ u^,{r') (59) 

Using the generic scaling law for S we find 

a;~rQ'';f~rQ S Tq " (60) 

We see, that the energy dissipation rate would stay finite in the limit of the vanishing 
viscous scale only for the Kolmogorov value of the index. This argument looks rather 
cheap, but I think it is basically correct. The constant value of the energy dissipation 
rate in the limit of vanishing viscosity arises as the quantum anomaly in the field theory, 
through the finite limit of the point splitting in the correspondent energy current .0 

There is another version of this argument, which I like better. The dynamics of Euler 
fluid in infinite system would not exist, for the non- Kolmogorov scaling. The extra powers 
of loop size would have to enter with the size L of the whole system, like (^) • So, in 
the regime with finite energy pumping rate S the infinite Euler system can exist only for 
the Kolmogorov index. This must be the essence of the original Kolmogorov reasoning 
i- 

The problem is that nobody proved that such limit exists, though. Within the usual 
framework, based on the velocity correlation functions, one has to prove, that the infrared 
divergencies, caused by the sweep, all cancel for the observables. Within our framework 
these problems disappear, as we shall see later. 

As for the correlation functions in inertial range, unfortunately those cannot be com- 
puted in the WKB approximation, since they involve the contour shrinking to a double 
line, with vanishing area inside. Still, most of the physics can be understood in loop terms, 
without these correlation functions. The large scale behavior of the loop functional reflects 
the statistics of the large vorticity structures, encircled by the loop. 

''I am grateful to A. Polyakov and E. Siggia for inspiring comments on this subject. 
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5 Loop Equation for the Circulation PDF 



The loop field could serve as the generating function for the PDF Pc{^) for the circulation. 
The Fourier integral 

dg 



c 



(61) 



can be analyzed in the same way as the loop field before. The only difference is that the 
factors of g appear in front of various terms. These factors can be replaced by 

d 



q ^ I 

^ dV 



(62) 



acting on -Pc(r)- 

As a result we find 
d_ 



Pc{T) = - 



d f 
dTJc 



dr„d, 



dVa 

C 

5Pc{T) 



-U[C] 



47r|r — r'l^ Sai3a(f)5ai3j(r') 



(63) 



All the imaginary units disappear, as they should. As for the viscosity and forcing, these 
terms can be neglected in inertial range in the same way as before. The only new thing 
is that one has to assume that F ^ u in inertial range in addition to above assumptions 
about the size of the loop. 

In absence of these terms there are no dimensional parameters so that the following 
scaling laws hold (with the same index k as before) 

^c(r) = t^'^-^F[^,^] (64) 

The factor t^'^"^ came from the normahzation of probability density. Note, that this is 
more general law than before. Here we do not have to use the WKB approximation for 
the PDF. In other words, the whole PDF rather than just its decay at large F satisfies 
the scaling law. 

The steady distribution would have the form of 



C 



where the scaling functional $ satisfies the homogeneous equation 



dr^ 



c 



dh' ''^ 



47r|r — r'p (5(7^a(r)5(T^^(r') 



= 



with the normalization condition 



r°° dr 

J-ooY 



- c 

K 



(65) 



(66) 



(67) 
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In principle, there could be different scaling functions for positive and negative F, 
rather than just absolute value |r| prescription. This would correspond to above men- 
tioned violation of the time reversal symmetry. However, as we mentioned above, there 
is no exact relation which would eliminate the symmetric solution. 

The Kolmogorov triple correlation function vanishes for vorticities (see Appendix I), so 
that there is no restriction on the asymmetric part of the circulation PDF. Nevertheless, 
the Kolmogorov scaling k = | seems to me the most likely possibility, by the reasons 
discussed in the previous section. 

The homogeneous loop equation requires some boundary conditions at large loops, to 
provide a meaningful solution. The asymptotic decrease of PDF 



Pc{r) ~ exp (-Q 
would lead to the same WKB equation as before 



Q^oo (68) 



.3./ r;-r^ SQ[C] SQ[C] 



d'r'-f L^_li^_^ = (69) 



We are studying this equation in the next section. 



6 Tensor Area law 

The Wilson loop in QCD decreases as exponential of the minimal area, encircled by the 
loop, leading to the quark confinement. What is the similar asymptotic law in turbulence? 
The physical mechanisms leading to the area law in QCD are absent here. Moreover, there 
is no guarantee, that ^[C] always decreases with the size of the loop. 

This makes it possible to look for the simple Anzatz, which was not acceptable in 
QCD, namely 

S[C] = s (Sj,) (70) 

where 

= / r^u (71) 
J c 

is the tensor area encircled by the loop C. The difference between this area and the 
scalar area is the positivity property. The scalar area vanishes only for the loop which 
can be contracted to a point by removal of all the backtracking. As for the tensor area, 
it vanishes, for example, for the 8 shaped loop, with opposite orientation of petals. 

Thus, there are some large contours with vanishing tensor area, for which there would 
be no decrease of the \1/ functional. In QCD the Wilson loops must always decrease at 
large distances, due to the finite mass gap. Here, the large scale correlations are known to 
exist, and play the central role in the turbulent flow. So, I see no convincing arguments 
to reject the tensor area Anzatz. 
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This Anzatz in QCD not only was unphysical, it failed to reproduce the correct short- 
distance singularities in the loop equation. In turbulence, there are no such singularities. 
Instead, there are the large-distance singularities, which all should cancel in the loop 
equation. 

It turns out, that for this Anzatz the (turbulent limit of the) loop equation is satisfied 
automatically, without any further restrictions. Let us verify this important property. 
The first area derivative yields 

/ \ ^ ds 

'^-M -J^)-'W. '''' 

The factor of 2 comes from the second term in the variation 

T TT — ^an^lSu — ^ai/^/J/i (73) 

Note, that the right side does not depend on r. Moreover, you can shift r aside from the 
base loop C, with proper wires inserted. The area derivative would not change, as the 
contribution of wires drops. 

This implies, that the corresponding vorticity uj^iy{r) is space independent, it only 
depends upon the loop itself. The velocity can be reconstructed from vorticity up to 
irrelevant constant sterms 



vpir) = -r^u';^f, (74) 



2 

This can be formally obtained from the above integral representation 



r„ — r' 



as a residue from the infinite sphere R= \r'\ oo. One may insert the regularizing factor 
\r'\~^ in a;, compute the convolution integral in Fourier space and check that in the limit 
e — > 0"*" the above linear term arises. So, one can use the above form of the loop equation, 
with the analytic regularization prescription. 

Now, the vuj term in the loop equation reads 

£ dr^ v'^ir) cx S^, c^g, u;^^ (76) 

This tensor trace vanishes, because the first tensor is antisymmetric, and the product of 
the last two antisymmetric tensors is symmetric with respect to 0:7. 

So, the positive and negative terms cancel each other in our loop equation, like the 
"income" and "outcome" terms in the usual kinetic equation. We see, that there is an 
equilibrium in our loop space kinetics. 

Prom the point of view of the notorious infrared divergencies in turbulence, the above 
calculation explicitly demonstrates how they cancel. By naive dimensional counting these 
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terms were linearly divergent. The space isotropy lowered this to logarithmic divergency 
in ([75|), which reduced to finite terms at closer inspection. Then, the explicit form of 
these terms was such, that they all cancelled. 

This cancellation originates from the angular momentum conservation in fiuid mechan- 
ics. The large loop C creates the macroscopic eddy with constant vorticity and linear 
velocity v^{r) oc r. This is a well known static solution of the Navier-Stokes equation. 
The eddy is conserved due to the angular momentum conservation. The only nontrivial 
thing is the functional dependence of the eddy vorticity upon the shape and size of the 
loop C. This is a function of the tensor area S^^, rather than a general functional of the 
loop. 

Combining this Anzatz with the space isotropy and the Kolmogorov scaling law, we 
arrive at the tensor area law 

m^e^v(^-B (^(S^;^)')') (77) 

The universal constant B here must be real, in virtue of the loop orientation symmetry. 
When the orientation is reversed C{9) C{2'k — 9), the loop integral changes sign, but 
its square, which enters here, stays invariant. Therefore, the constant in front must be 
real. The time reversal tells the same, since both viscosity u and the energy dissipation 
rate £ are time-odd. Therefore, the ratio ^ is time-even, hence it must enter \1/[C] with 
the real coefficient. Clearly, this coefficient B must be positive, since |\I'[C]| < 1. 

Note, however, that we did not prove this law. The absence of decay for large twisted 
loops with zero tensor area is suspicious. Also, the physics seems to be different from 
what we expect in turbulence. The uniform vorticity, even a random one, as in this 
solution, contrasts the observed intermittent distribution. Besides, there clearly must be 
corrections to the asymptotic law, whereas the tensor area law is exact. This is far too 
simple. We discussed this unphysical solution mostly as a test of the loop technology. 



7 Scalar Area law 

Let us now study the scalar area law, which is a valid Anzatz for the asymptotic decay 
of the circulation PDF. The set of equations for the minimal surface is summarized in 
Appendix A. All we need here is the following representation 

l/"/" fix — v)'^\ 

^^2lJ J J '^^^^^(^)'^^^^^'(y) — Jl — I (78) 

where Lr = \T\i£^^. The distance (x — yY is measured in 3-space and integration goes 
along the minimal surface. It is implied that its size is much larger than L^. 

In this limit the integration over, say, y can be performed along the local tangent plane 
at X in small vicinity y — x ^ L-p , after which the factors of cancel. We are left then 
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with the ordinary scalar area integral 



^ J da^^{x)da^^{y)6^{x - v) ^ J d^^y/g (79) 



A 

2 

In the previous, regularized form the area represents so called Stokes functional^ 
which can be substituted into the loop equation. The area derivative of the area reads 



In the local limit this reduces to the tangent tensor 

5A 



da^^{y)6 {x - y) = t^^{x) (81) 



It is implied that the point x approaches the contour from inside the surface, so that the 
tangent tensor is well defined 

Here is the local tangent vector of the loop, and riy is the inside normal to the loop 
along the surface. 

The second area derivative of the regularized area in this limit is just the exponential 



5^ A _ 1 / {x-yf 
6aap{x)5a^5{y) Ll \ ^ Ll 



(83) 



Should we look for the higher terms of the asymptotic expansion at large area we would 
have to take into account the shape of the minimal surface, but in the thermodynamical 
limit we could neglect the curvature of the loop and use the planar disk. 

Let us use the general WKB form of PDF 

^c(r) = ^exp(-g(A,_i_)) (84) 

We shall skip the arguments of effective action Q. We find on the left side of the loop 
equation 

OtQdrQ - dtdrQ (85) 
On the right side we find the following integrand 

{{d^Qf - dlQ) - dAQ , (86) 

The last term drops after the r' integration in virtue of symmetry. The leading terms 
in the WKB approximation on both sides are those with the first derivatives. We find 

OtQdrQ = {OaQ? / dr^-^ f d'r' ~ % (87) 
Jc 0aai3[r) J 47r|r — r'l'^ da^fs[r) 

18 



In the last integral we substitute above explicit form of the area derivatives and perform 
the d^r' integration first. In the thermodynamical limit only the small vicinity r' — y '-^ Ly 
contributes, and we find 



J At^W - r'\^ Sa^nir') ^^^^^ 



47r|r — r'\^ Sa^p{r') J 47r|r — 

This integral logarithmically diverges. We compute it with the logarithmic accuracy 
with the following result 



1. ^ 



oc ^ In (89) 



47r|r — y\^ tt A 

The meaning of this integral is the average velocity in the WKB approximation. This 
velocity is tangent to the loop, up to the next correction terms at large area. 

Now, the emerging loop integral vanishes due to symmetry 



^ dratptaf: 



(90) 



as the line element dr^ is directed along the tangent vector ta, and the tangent tensor 
tap is antisymmetric. Similar mechanism was used in the tensor area solution, only there 
the cancellations emerged at the global level, after the closed loop integration. Here the 
right side of the loop equation vanishes locally, at every point of the loop. Anyway, we 
see, that the scalar area indeed represents the steady solution of the loop equation in the 
leading WKB approximation. 

It might be instructive to compare this solution with another known exact solution of 
the Euler dynamics, namely the Gibbs solution 

P[v] = exp J d'r^vl^ (91) 

For the loop functional it reads 

*c(7) = exp (-^ £ dra £ dr'f^S^r - r')) (92) 

The integral diverges, and it corresponds to the perimeter law 

j) dva <j> dr'gS^{r — r') Tq ^ j> \dr\ (93) 
where tq is a small distance cutoff. For the PDF it yields 

When the Gibbs solution is substituted into the loop equation, we observe the same 
thing. Average velocity is tangent to the loop, which leads to vanishing integrand in the 
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loop equation. The difference is tliat in our case tliis is true only asymptotically, there 
are next corrections. 

The shape of the function Q is not fixed by this equation in the leading WKB approx- 
imation. In a scale invariant theory it is natural to expect the power law 

^(jl.isJiT)-"-' (r"^'-")" (95) 

There is one more arbitrary index /i involved. Even for the Kolmogorov law = | the F 
dependence remains unknown. 



8 Discussion 

So, we found two asymptotic solutions of the loop equation in the turbulent limit, not 
counting the Gibbs solution. It remains to be seen, which one (if any) is realized in 
turbulent flows. The tensor area solution is mathematically cleaner, but its physical 
meaning contradicts the intermittency paradigm. It corresponds to the uniform vorticity 
with random magnitude and random direction, rather that the regions of high vorticity 
interlaced with regions of low vorticity, observed in the turbulent flows. 

The recent numerical experiments^] favor the scalar area rather than the tensor one. 
Also, the Kolmogorov scaling was observed in these experiments. The Reynolds number 
was only ~ 100 which was too small to make any conclusions. We have to wait for the 
experiments ( real or numerical ) with the Reynolds numbers few orders of magnitude 
larger. 

The scalar area is less trivial than the tensor one. The minimal functional of 

the loop cannot be represented as any explicit contour integral of the Stokes type, therefore 
it corresponds to infinite number of higher correlation functions present. Moreover, there 
could be several minimal surfaces for the same loop, as the equations for the minimal 
surface are nonlinear. Clearly, the one with the least area should be taken. 

The natural generalization of this solution is the string Anzatz where the sum over all 
surfaces bounded by the loop is taken 

Pci^)= E -p(-g(A,_L_)) (96) 

S'.dS — C 

At large loop the minimal Q terms will remain. The extremum condition 

5Q = ^5 A = (97) 

will be satisfied for the minimal surface. 

However, the sum over random surfaces is not well defined. The recent studies]^ 
indicate that the typical closed surfaces degenerate to branched polymers. For the surface 
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bounded by a fixed loop this cannot happen, of course. Still nobody knows how to 
compute such sums. The loop equation in principle allows to systematically compute the 
corrections to the area law as the WKB expansion. 

The WKB solution is incomplete so far. The leading term in the loop equation is 
annihilated by arbitrary function of the area (scalar or tensor). The similar ambiguity 
was present in the Gibbs solution, where arbitrary function of the Hamiltonian satisfied 
the Liouville equation for the velocity PDF. In that case the ambiguity was removed by 
extra requirement of thermodynamic limit: only the exponential of the hamiltonian would 
agree with the factorization of the PDF for two remote parts of the system. 

What could be a similar requirement here? The area of the minimal surface represents 
the effective volume of the system at large loop. The circulation can be written as a surface 
integral of vorticity, which makes the circulation an extensive variable at this surface. The 
average vorticity 

r . . 

= - 98 

A ^ ' 

represents an intensive quantity. The thermodynamic limit would then correspond to 

Q = Aq{u) (99) 

Comparing this with the previous formula for Q (|95|) we conclude that /i = 1. In this 

case 

Q = const Aoo'^'' (100) 

In principle, there could be two different laws for positive and negative P, due to violation 
of the time reversal invariance 

Q q±A\cu\'^'' (101) 

Another line of argument might start with an assumption of decorrelated average 
vorticity Ui at various parts Si of the area Aq of the minimal surface. The net circulation, 
adding up from the large number n ~ ^ 1 of independent random terms UiAq would 
be a gaussian variable as a consequence of the law of large numbers. We would have then 



Q JT^ = -nr^ (102) 

nAlujf AAouf 



This would agree with the previous estimate at 

1 



so that 



H = -,Aoiof A^'^- (103) 

K 



Q const T'^A^-^ (104) 



The natural assumption here would be that the vorticity variance Ui does not scale with 
the area A, so that 

Aq ~ A^-^ (105) 
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The Gaussian behavior (with k 



I ) was observed in numerical experiments 0, but 



the Reynolds number was too low to make conclusions at this point. There could be a 
scahng function 



which starts linearly and then grows as a power, say, q{x) = (1 + ax)'^. I suggest that 
this function should be studied in real and numerical experiments. This would teach us 
something new about turbulence. 
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A Loop Expansion 



Let us outline the method of direct iterations of the loop equation. The full description 
of the method can be found in The basic idea is to use the following representation 
of the loop functional 

nC] = l + f:Uidr^K..ldr:"] W:^,,,^Jn,...rn) (107) 
n=2 ^ Jc ) cychc 

This representation is valid for every translation invariant functional with finite area 
derivatives (so called Stokes type functional). The coefficient functions W can be related 
to these area derivatives. The normalization \E'[0] = 1 for the shrunk loop is implied. 

In general case the integration points ri, . . .r„ in ( |107|) are cyclicly ordered around 
the loop C . The coefficient functions can be assumed cyclicly symmetric without loss 
of generality. However, in case of fluid dynamics, we are dealing with so called abelian 
Stokes functional. These functionals are characterized by completely symmetric coefficient 
functions, in which case the ordering of points can be removed, at expense of the extra 
symmetry factor in denominator 

nc] = i + Y.hi ■■■i <"^<5^i..-. (^1' • • • (108) 

n=2 ■"^ •'^ 

The incompressibility conditions 

5a.W/:^...a„(n,...r„) = (109) 

does not impose any further restrictions, because of the gauge invariance of the loop 
functionals. This invariance (nothing to do with the symmetry of dynamical equations!) 
follows from the fact, that the closed loop integral of any total derivative vanishes. So, 
the coefficient functions are defined modulo such derivative terms. In effect this means, 
that one may relax the incompressibility constraints ( |109| ), without changing the loop 
functional. 

To avoid confusion, let us note, that the physical incompressibility constrains are not 
neglected. They are, in fact, present in the loop equation, where we used the integral 
representation for the velocity in terms of vorticity. Still, the longitudinal parts of W 
drop in the loop integrals. 

The loop calculus for the abelian Stokes functional is especially simple. The area 
derivative corresponds to removal of one loop integration, and differentiation of the cor- 
responding coefficient function 

y^J^lZ...Ar,r^,...T^) (110) 
d.ho^ (111) 



^<yiiv{r) ^^n\Jc "Jc 
where 
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In the nonabelian case, there would also be the contact terms, with W at coinciding 
points, coming from the cyclic ordering 0. In abelian case these terms are absent, since 
W is completely symmetric. 

As a next step, let us compute the local kinetic term 

(112) 



Jc ocT^u[r) 



Using above formula for the loop derivative, we find 



(r, ri, . . .r„) 



(113) 



The net result is the second derivative of W with respect to one variable. Note, that the 
second term in dropped, as the total derivative in the closed loop integral. 

As for the nonlocal kinetic term, it involves the second area derivative off the loop, 
at the point r', integrated over r' with the corresponding Green's function. Each area 
derivative involves the same operator V , acting on the coefficient function. Again, the 
abelian Stokes functional simplifies the general framework of the loop calculus. The 
contribution of the wires cancels here, and the ordering does not matter, so that 



OO 1 /• r I 



(114) 



Using these relations, we can write the steady state loop equation as follows 

"9. 



Lb- 



in+V 



+ 





Fig. 6 



Here the light dotted lines symbolize the arguments ak, of W , the big circle denotes 
the loop C, the tiny circles stand for the loop derivatives, and the pair of lines with 
the arrow denote the Green's function. The sum over the tensor indexes and the loop 
integrations over are implied. 

The first term is the local kinetic term, the second one is the nonlocal kinetic term, 
and the right side is the potential term in the loop equation. The heavy dotted line in 
this term stands for the correlation function F of the random forces. Note that this term 
is an abelian Stokes functional as well. 
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The iterations go in the potential term, starting with ^[C] = 1. In the next approx- 
imation, only the two loop correction W^^^,^{ri,r2) is present. Comparing the terms, we 
note, that nonlocal kinetic term reduces to the total derivatives due to the space symme- 
try (in the usual terms it would be {vu>) at coinciding arguments), so we are left with the 
local one. 



This yields the equation 

u'd''W^f,ir-r')^F{r-r')5a/3 
modulo derivative terms. The solution is trivial in Fourier space 



r') 



d^k F{k) 
^exp(zMr-r))5«,^ 



Note, that we did not use the transverse tensor 



kakp 



(115) 



(116) 



(117) 



Though such tensor is present in the physical velocity correlation, here we may use Sa/s 
instead, as the longitudinal terms drop in the loop integral. This is analogous to the 
Feynman gauge in QED. The correct correlator corresponds to the Landau gauge. 

The potential term generates the four point correlation FW'^. which agrees with the 
disconnected term in the on the left side 



(ri,r2,r3,r4) 



^2) Wl^^ (rs - + 



(118) 



<«3 (^1 - ^3) W^4a4 (^2 - n) + wl^^ in - n) wl^^,^^ {r, - r,) 



47r|) — r'p ■ 



In the same order of the loop expansion, the three point function will show up. The 
corresponding terms in kinetic part must cancel among themselves, as the potential term 
does not contribute. The local kinetic term yields the loop integrals of d'^W^, whereas the 
nonlocal one yields VW"^ V'W"^, integrated over d'^r' with the Greens's function ^ 
The equation has the structure 



= 





Fig. 7 



Now it is clear, that the solution of this equation for would be the same three point 
correlator, which one could obtain (much easier!) by direct iterations of the Navier-Stokes 
equation. 
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The purpose of this painful exercise was not to give one more method of developing 
the expansion in powers of the random force. We rather verified that the loop equations 
are capable of producing the same results, as the ordinary chain of the equations for the 
correlation functions. 

In above arguments, it was important, that the loop functional belonged to the class 
of the abelian Stokes functionals. Let us check that our tensor area Anzatz 



yC 



c 



Tadr 



belongs to the same class. Taking the square we find 



c Jc ^ 



1 



c 



dr'Jr — r')^ 



c 



(119) 



(120) 



where the last transformation follows from the fact, that only the cross term in (r — r'Y 
yields nonzero after double loop integration. 

Any expansion in terms of the square of the tensor area reduces, therefore to the 
superposition of multiple loop integral of the product of (r^ — Tj)^, which is an example of 
the abelian Stokes functional. In the limit of large area, this could reduce to the fractional 
power. An example could be, say 





( 




^[C] = exp 


B 


H 









(121) 



One could explicitly verify all the properties of the abelian Stokes functional. This example 
is not realistic, though, as it does not have the odd terms of expansion. In the real world 
such terms are present at the viscous scales. According to our solution, this asymmetry 
disappears in inertial range of loops (which does not apply to velocity correlators at 
inertial range, as those correspond to shrunk loops). 



B Matrix Model 

The Navier-Stokes equation represents a very special case of nonlinear PDE. There is a 
well known galilean invariance 

Va{r,t) ^ Va{r - ut,t) + Ua (122) 

which relates the magnitude of velocity field with the scales of time and space. Let us 
make this relation more explicit. 

^At the same time it tells us that the constant part of velocity if frame dependent, so that it better 
be eliminated, if we would like to have a smooth limit at large times. Most of notorious large scale 
divergencies in turbulence are due to this unphysical constant part. 
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First, let us introduce the vorticity field 

i^nu = d^Vu - d^v^ (123) 
and rewrite the Navier-Stokes equation as follows 

Va = i^dpupa - vpupa " d^w; w^p+— (124) 



This w is the well known enthalpy density, to be found from the incompressibility 
condition divv = 0, i.e. 

d^W = daVfsUfSa (125) 

As a next step, let us introduce "covariant derivative" operator 

Da = l^da - (126) 

and observe that 

2 [DaDp] = i^cupa (127) 

2Dp [DaDp] + h.c.^ T^dpUfSa - vpUpa (128) 

where h.c. stands for hermitean conjugate. 

These identities allow us to write down the following dynamical equation for the 
covariant derivative operator 

Da = Dp [DaDp] - DaW + h.c. (129) 

As for the incompressibility condition, it can be written as follows 

DaDi] = (130) 
The enthalpy operator = ^ is to be determined from this condition , or, equivalently 

[Da [DaW]] = [Da, Dp [DaDp]] (131) 

We see, that the viscosity disappeared from these equations. This paradox is resolved 
by extra degeneracy of this dynamics: the antihermitean part of the D operator is con- 
served. Its value at initial time is proportional to viscosity. 

The operator equations are invariant with respect to the time independent unitary 
transformations 

Da S'fDaS; S^S = 1 (132) 
and, in addition, to the time dependent unitary transformations with 

S{t) = exp (^^tup (Dp - Dt)) (133) 
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corresponding to the galilean transformations. 
One could view the operator as the matrix 



(134) 



where the functions ipj{r) are the Fourier of Tchebyshev functions depending upon the 
geometry of the problem. 

The finite mode approximation would correspond to truncation of this infinite size 
matrix to finite size A^. This is not quite the same as leaving terms in the mode 
expansion of velocity field. The number of independent parameters here is 0{N'^) rather 
then 0{N). It is not clear whether the unitary symmetry is worth paying such a high 
price in numerical simulations! 

The matrix model of Navier-Stokes equation has some theoretical beauty and raises 
hopes of simple asymptotic probability distribution. The ensemble of random hermitean 
matrices was recently applied to the problem of Quantum Gravity 0, which led to a 
genuine breakthrough in the field. 

Unfortunately, the model of several coupled random matrices, which is the case here, 
is much more complicated then the one matrix model studied in Quantum Gravity. The 
dynamics of the eigenvalues is coupled to the dynamics of the "angular" variables, i.e. 
the unitary matrices S in above relations. We could not directly apply the technique of 
orthogonal polynomials, which was so successful in the one matrix problem. 

Another technique, which proved to be successful in QCD and Quantum Gravity is 
the loop equations. This method, which we are discussing at length in this paper, works 
in field theory problems with hidden geometric meaning. The turbulence proves to be an 
ideal case, much simpler then QCD or Quantum Gravity. 

C The Reduced Dynamics 

Let us now try to reproduce the dynamics of the loop field by a simpler Anzatz 



The difference with original definition (|^) is that our new function Pa{0) depends directly 
on 6 rather then through the function t^a(r) taken at Tq, = Ca{d). This is the d ^ 1 
dimensional reduction we mentioned before. From the point of view of the loop functional 
there is no need to deal with field v{r) , one could take a shortcut. 

Clearly, the reduced dynamics must be fitted to the Navier-Stokes dynamics of original 
field. With the loop calculus, developed above, we have all the necessary tools to build 
this reduced dynamics. 




(135) 
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Let us assume some unknown dynamics for the P field 



p^{e) = F^{e,[p]) (136) 

and compare the time derivatives of original and reduced Anzatz. We find in (|135|) instead 
of (H 

'-jdCMFAOAP]) (137) 

Now we observe, that P' could be replaced by the functional derivative, acting on the 
exponential in (|135|) as follows 



-tvP'^{e) (138) 



This means, that one could take the operators of the Section 2, expressing velocity and 
vorticity in terms of the spike operator, and replace the functional derivative as above. 
This yields the following formula for the spike derivative 

D^{Q, t) = -ivj^ # h - ' ^' j P'M = ^ \ ^ sgn(/i)P„ {e + e(l + ^)) 

(139) 

This is the weighted discontinuity of the function P{9), which in the naive limit e — 
would become the true discontinuity. However, the function P{9) has in general the 
stronger singularities, then discontinuity, so that this limit cannot be taken yet. 

Anyway, we arrive at the dynamical equation for the P field 

= uDp^lp^ - Vp^p^ (140) 

where the operators V, VL of the Section 2 should be regarded as the ordinary numbers, 
with definition (|139| ) of D in terms of P. 

All the functional derivatives are gone! We needed them only to prove equivalence of 
reduced dynamics to the Navier-Stokes dynamics. 

The function Pa{0) would become complex now, as the right side of the reduced 
dynamical equation is complex for real Pa{6)- 

Let us discuss this puzzling issue in more detail. The origin of imaginary units was 
the factor of z in exponential of the definition of the loop field. We had to insert this 
factor to make the loop field decreasing at large loops as a result of oscillations of the 
phase factors. Later this factor propagated to the definition of the P field. 

Our spike derivative D is purely imaginary for real P, and so is our Q operator. This 
makes the velocity operator V real. Therefore the DQ term in the reduced equation (|140|) 
is real for real P whereas the VQ term is purely imaginary. 

This does not contradict the moments equations, as we saw before. The terms with 
even/odd number of velocity fields in the loop functional are real/imaginary, but the 
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moments are real, as they should be. The complex dynamics of P simply doubles the 
number of independent variables. 

There is one serious problem, though. Inverting the spike operator we implicitly 
assumed, that it was antihermitean, and could be regularized by adding infinitesimal 
negative constant to D"^ in denominator. This, indeed, works perturbatively, in each term 
of expansion in time, or that in size of the loop, as we checked. However, beyond this 
expansion there would be a problem of singularities, which arise when -0^(6*) vanishes at 
some 9. 

In general, this would occur for complex 9, when the imaginary and real part of D^{9) 
simultaneously vanish. One could introduce the complex variable 

e^&^z; e-'' = -; cf d9 = f - (141) 

where the contour of z integration encircles the origin around the unit circle. Later, in 
course of time evolution, these contours must be deformed, to avoid complex roots of 
DUO)- 



D Initial Data 

Let us study the relation between the initial data for the original and reduced dynamics. 
Let us assume, that initial field is distributed according to some translation invariant 
probability distribution, so that initial value of the loop field does not depend on the 
constant part of C{9). 

One can expand translation invariant loop field in functional Fourier transform 

*[C] = j DQ5'(^^d4>Q{4>)^W[Q]exp(^i^d9C^{9)Q^{9)^ (142) 

which can be inverted as follows 

(^j> d(i>Q{^)^ W[Q] = j DC^![C\ exp [-i j> d9C^{9)Q^{e)^ (143) 

Let us take a closer look at these formal transformations. The functional measure for 
these integrations is defined according to the scalar product 

r rlf) 

iA,B)^f-A{9)B{9) (144) 
which diagonalizes in the Fourier representation 

+ 00 

Ai9) = J2Ane"''; A_^ = Al (145) 



+ 00 



{A, B)^Y. ^nB-n = AoBo + ^ a'X + <K; < = ^^K. < = v^^^n (146) 
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The corresponding measure is given by an infinite product of the Euchdean measures 
for the imaginary and real parts of each Fourier component 

oo 

DQ = d'QoUd'q'nd'q'n (147) 
1 

The orthogonahty of Fourier transformation could now be explicitly checked, as 

J DCew (^ J deC^ie) {A^{e) - B^{e))^ (148) 

= J d'Co n d'c'jX exp (27rz {A^ - B,) +± < « - b'J + < « - b:)^ ) 

= SHAo-Bo)l[^'{a:,-b'J6^a:-b:) 
1 

Let us now check the parametric invariance 

e ^ f{ey, f{2n) - /(O) = 2n; f{9) > (149) 

The functions C{9) and P{9) have zero dimension in a sense, that only their argument 
transforms 

c{e)^c{f{e))- p{e)^p{f{9)) (iso) 

The functions Q{9) and P'{0) in above transformation have dimension one 

p'{e) ^ f{e)p' im) ; Q{e) ^ f{e)Q (m) (isi) 

so that the constraint on Q remains invariant 

jdeQ{e) = jdf{e)Q{f{e)) (152) 

The invariance of the measure is easy to check for infinitesimal reparametrization 

/(6») + e(6'); e(27r) - e(0) (153) 

which changes C and (C, C) as follows 

r df) f df) 

5C{9) = e{e)C\9); S{C, C) ^ j -e{e)2C^{9)CM ^ - j> -e\9)Cl{e) (154) 

The corresponding Jacobian reduces to 

l-jd0e\9)^l (155) 

in virtue of periodicity. 

This proves the parametric invariance of the functional Fourier transformations. Using 
these transformations we could find the probability distribution for the initial data of 

Pa{e)^-y f d(l>QM (156) 
Jo 
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The simplest but still meaningful distribution of initial velocity field is the Gaussian 
one, with energy concentrated in the macroscopic motions. The corresponding loop field 
reads 

vi>o[q = exp(^--^dc„(^) jdc^{e')f{c{e)-c{e'))^ (157) 

where f{r — r') is the velocity correlation function 

{vair)vpir')) = [6^^ - d^O^d;') /(r - r') (158) 

The potential part drops out in the closed loop integral. 

The correlation function varies at macroscopic scale, which means that we could ex- 
pand it in Taylor series 

/(r-r')^/o-/i(r-r')' + ... (159) 
The first term /q is proportional to initial energy density, 

I {vl) = ^/o (160) 
and the second one is proportional to initial energy dissipation rate 

So = -V {v^dlv^) = 2d{d - l)z./i (161) 

where d = 3 is dimension of space. 

The constant term in ( |159| ) as well as + r'^ terms drop from the closed loop integral, 
so we are left with the cross term rr' 

^o[c] ^ exp (^-/i fdC49) j>dc^{e')Cp{e)Cp{e')^ (162) 

This is almost Gaussian distribution: it reduces to Gaussian one by extra integration 

^o[C] ^ const j c/^cuexp (-ul^) exp {2i^^u^, <f dC^(^)a(^)) (163) 

The integration here goes over all ^i^tli = 3 independent a < P components of the 
antisymmetric tensor a;^^. Note, that this is ordinary integration, not the functional one. 
The physical meaning of this u is the random constant vorticity at initial moment. 

At fixed uj the Gaussian functional integration over C 

J DCexp [ijde {^Cp{e)p'(,{e) + 2/a^„/3c;(^)c^(^))) (i64) 

can be performed explicitly, it reduces to solution of the saddle point equation 

P'p{e) = AuJhup^C'M (165) 
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which is trivial for constant uj 

Ca{e) = ^^a;-jP^(^) (166) 

The inverse matrix is not* unique in odd dimensions, since Det u>ap — 0. However, the 
resulting pdf for P is unique. This is the Gaussian probability distribution with the 
correlator 

{Pa{9)Pf3{9')) = 2iv^^u^psign{e' - 9) (167) 

Note, that antisymmetry of u compensates that of the sign function, so that this 
correlation function is symmetric, as it should be. However, it is antihermitean, which 
corresponds to purely imaginary eigenvalues. The corresponding realization of the P 
functions is complex! 

Let us study this phenomenon for the Fourier components. Differentiating the last 
equation with respect to 6 and Fourier transforming we find 




(168) 



This cannot be realized at complex conjugate Fourier components P„,-n = Pan t*ut 
we could take Pa,n = Pa,-n and Pa,n as real random variables, with correlation function 

{Pa,nP0,m) = ^ ^nm ; n > (169) 

The trivial realization is 

4z/ / — 

Pa,n = —\J h(^al3Pp,n (170) 

with Pp^n being Gaussian random numbers with unit dispersion. 

As for the constant part P^.o of Pa{9) , it is not defined, but it drops from equations 
in virtue of translational invariance. 



E W-functional 

The difficulties of turbulence are hidden in the loop equation, but they show up, if you 
try to solve it numerically. The main problem is that one cannot get rid of the cutoffs 
e, (5 in the definitions of the spike derivatives. These cutoffs are designed to pick up 
the singular contributions in the angular integrals, but with finite number of modes, such 
as Fourier harmonics there would be no singularities. We did not find any way to truncate 
degrees of freedom in the P equation, without violating the parametric invariance. It very 
well may be, that this invariance would be restored in the limit of large number of modes, 
but it looks that there are too much ambiguity in the finite mode approximation. 
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After some attempts, we found the simpler version of the loop functional, which can 
be studied analytically in the turbulent region. This is the generating functional for the 
scalar products Pa{(^i) Pa{(^2) 



W[S] = ^exp f dOi f de2S{9i,92)Pa{ei)Pa{92: 

where, as before, the averaging goes over initial data for the P field. 
The time derivative of this W-functional 

w = -2l<lde, <fde2S{9^,e2)PMPa{02)exp (- Ue, i de2S{d^,d2)PMP^{e2 



[in) 



can be expressed in terms of functional derivatives of W by replacing 

5 



Pai4>l)Pc 



aW2 



6S{(f)i, (p2 



(172) 

(173) 



for every scalar product of P fields, which arise after expansion of the spike derivatives 
( |139| ), ([TE|), ( pup in the scalar product 

Pa{ei)Pa{02) = I^Pa{0l)Df,{e2)npa{02) - P. (^1 ) l^M^2)^^/3a (^2) (174) 

This equation has the structure 

w= id^esiei,e2) (a^ 



' 5 ' 




' 5 ' 




_6S_ 


W + A3 


6S_ 





(175) 

where A^. [X] stands for the k— degree homogenous functional of the function X{6i, 62). 

The operator is also the homogeneous functional of the negative degree k = —1. 
It can be written as follows 



with 



POO 

D-\e.e)W\S]= drWlS + rU] 
Jo 

U {61,62) = e-^ sgn{e + e-6i) sgn(6' + 6-62] 



(176) 
(177) 



F Possible Numerical Implementation 

The above general scheme is fairly abstract and complicated. Could it lead to any practical 
computation method? This would depend upon the success of the discrete approximations 
of the singular equations of reduced dynamics. 

The most obvious approximation would be the truncation of Fourier expansion at 
some large number A^. With Fourier components decreasing only as powers of n this 
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approximation is doubtful. In addition, such truncation violates the parametric invariance 
which looks dangerous. 

It seems safer to approximate P{0) by a sum of step functions, so that it is piecewise 
constant. The parametric transformations vary the lengths of intervals of constant P{0), 
but leave invariant these constant values. The corresponding representation reads 

N 

PaiO) = E M + 1) - Po.il)) Q{9- 6,) ■ p{N + 1) = p{l), p{0) = (178) 

1=0 

It is implied that Oq = < 9i < 02 ■ ■ ■ < 0^ < 27r. By construction, the function P{6) 
takes value p{l) at the interval 9i^i < 9 < 9i. 

We could take P{9) at the middle of this interval as approximation to p{l). 

p{l)^P{9,)- e, = ]^{9i^, + 9{) (179) 
As for the time evolution of angles 9i , one could differentiate (|1 78|) in time and find 

Pa{9) = E iPM + 1) - l^-(O) Q{0-9i)-jz (P»(^ + 1) - P-(^)) - (180) 
;=o 1=0 

from which one could derive the following approximation 

n _ (Pa(0 -Pa(^ + 1)) [^'+' .nh (m f^s^^ 

- ( 7TTT2 /, d9P,,{9) (181) 

{p^{l + l) -p^il)) Jei 

The extra advantage of this approximation is its simplicity. All the integrals involved 
in the definition of the spike derivative ( |139| ) are trivial for the stepwise constant P{9). 
So, this approximation can be in principle implemented at the computer. This formidable 
task exceeds the scope of the present work, which we view as purely theoretical. 



G Uniqueness of the tensor area law 

Let us address the issue of the uniqueness of the tensor area solution. Let us take the 
following Anzatz 

S[C] = f l^j^dro,j>^dr'^W{r-r')^ (182) 

When substituted into the static loop equation (with the area derivatives computed in 
Appendix A), it yields the following equation for the correlation function W{r) 

= i dvo^i dr'^ i dr'!^Ua,p^{r, r', r") (183) 

J C " C " c 

UaiS'jir, r', r") = W{r - r')Vp^W{r' - r") + permutations 
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The derivative /' of the unknown function drops from the static equation. 

This equation should hold for arbitrary loop C. Using the Taylor expansion for the 
Stokes type functional 0, we can argue, that the coefficient function U must vanish up 
to the total derivatives. An equivalent statement is that the third area derivative of this 
functional must vanish. Using the loop calculus (see Appendix A) we find the following 
equation 

= Kuy'tu'y'W'Uc.a'a''{r, r', r") (184) 

which should hold for arbitrary r, r', r" . This leads to the overcomplete system of equations 
for W{r) in general case. However, for the special case W{r) = which corresponds to 
the square of the tensor area 

the system is satisfied as a consequence of certain symmetry. In this case we find in the 
loop equation 

2 j) dva j) dr'pijr — r'Y <j> dr'^ (^r^ — r^^ oc Sj^^ <j> dr^ j> dr'p{r — r'Y (186) 

The last integral is symmetric with respect to permutations of a, whereas the first 
factor is antisymmetric, hence the sum over a/5 yields zero, as we already saw above. 

It was assumed in above arguments, that the loop C consist of only one connected 
part. Let us now consider the more general situation, with arbitrary number n of loops 
Ci, . . . C„. The corresponding Anzatz would be 

S'„[Ci,...C„] = s„(si,...S") (187) 



where S* are tensor areas. 

This function should obey the same WKB loop equations in each variable. Introducing 
the loop vorticities 

(IBS) 

which are constant on each loop, we have to solve the following problem. What are the 
values of u;^^ such that the single velocity field v^ij-) could produce them? 

We do not see any other solutions, but the trivial one, with all equal a;^^, and linear 
velocity, as before. This would correspond to 

n „ 

s„ (S\ . . . S") = (S) ; S^. = E = f ^,dr, (189) 



k=l 



The loop equation would be satisfied like before, with C = WCfc. This corresponds to the 
additivity of loops 

Sr,[Cu...Cn] = S,[\£Ck] (190) 
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Note, that such additivity is the opposite to the statistical independence, which would 
imply that 

S^[C,,...Cn]=J2Si[Ck] (191) 

The additivity could also be understood as a statement, that any set of n loops is equiva- 
lent to a single loop for the abelian Stokes functional. Just connect these loops by wires, 
and note that the contribution of wires cancels. So, if the area law holds for arbitrary 
single loop, than it must be additive. 

This assumption may not be true, though, as it often happens in the WKB approxi- 
mation. There is no single asymptotic formula, but rather collections of different WKB 
regions, with quantum regions in between. In our case, this corresponds to the following 
situation. 

Take the large circular loop, for which the WKB approximation holds, and try to split 
it into two large circles. You will have to twist the loop like the infinity symbol oo, in 
which case it intersects itself. At this point, the WKB approximation might break, as 
the short distance velocity correlation might be important near the self-intersection point. 
This may explain the paradox of the vanishing tensor area for the oo shaped loop. Prom 
the point of view of our area law such loop is not large at all. 



H Minimal surfaces 



Let us present here the modern view at the classical theory of the minimal surfaces. The 
minimal surface can be described by parametric equation 

^:r« = X,(ei,6) (192) 
The function Xa{C) should provide the minimum to the area functional 

A[X] = ^rf^ = I d'^VD^ (193) 

where 

Gab = daX,d,X,, (194) 

is the induced metric. For the general studies it is sometimes convenient to introduce the 
unit tangent tensor as an independent field and minimize 

A [X, t,X] = J d^i ( eabdaX^d,X, ^ + A (l - tjj) (195) 

Prom the classical equations we will find then 

t,,= ^^daX^d,X,-tl, = l, (196) 

which shows equivalence to the old definition. 



37 



For the actual computation of the minimal area it is convenient to introduce the 
auxiliary internal metric gab 

A [X, g] = ll d^i tr g-^G ^D^. (197) 
The straightforward minimization with respect to gab yields 

gabiig-^G = 2Gab, (198) 

which has the family of solutions 

gab = XGab. (199) 

The local scale factor A drops from the area functional, and we recover original definition. 
So, we could first minimize the quadratic functional (|197|) with respect to X{^) (the linear 
problem), and then minimize with respect to gab (the nonlinear problem). 

The crucial observation is the possibility to choose conformal coordinates, with the 
diagonal metric tensor 

5, 



9ab = Sabp, 9ab = V ^ct g = p] (200) 
after which the local scale factor p drops from the integral 

A[X, P] = lj^ (f^daX^daX^. (201) 

However, the p field is implicitly present in the problem, through the boundary conditions. 

Namely, one has to allow an arbitrary parametrization of the boundary curve C. We 
shall use the upper half plane of f for our surface, so the boundary curve corresponds to 
the real axis ^2 = 0. The boundary condition will be 

X^(6,+0) = C(/(6)), (202) 

where the unknown function f{t) is related to the boundary value of p by the boundary 
condition for the metric 

gii = P = Gn = {d,X,f = C'^f" (203) 

As it follows from the initial formulation of the problem, one should now solve the linear 
problem for the X field, compute the area and minimize it as a functional of /(.). As we 
shall see below, the minimization condition coincides with the diagonality of the metric 
at the boundary 

[d,X,d2X,]^^^^^ = (204) 

The linear problem is nothing but the Laplace equation d^X = in the upper half 
plane with the Dirichlet boundary condition (|202D. The solution is well known 
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The area functional can be reduced to the boundary terms in virtue of the Laplace equa- 

d6 [X^d^X,]..^ (206) 



Substituting here the solution for X we find 



+00 r+oo 



dt 



dt 



,C ,{f{t))C,{f{t')) 

{t-t'-ioy 



This can be rewritten in a nonsingular form 



A[f] 



1 

An 



dt 



{t- t'f 



(207) 



(20^ 



which is manifestly positive. 

Another nice form can be obtained by integration by parts 

1 c+oo r+oo 

m = ^ / dtnt) / dt' nt')c'^{m) cum) log \t - 1'\ (209) 

ZTT J —00 J —00 

This form allows one to switch to the inverse function t(/) which is more convenient for 
optimization 



1 r+oo r+oo 

^[^^^2^ J-00 ^^/-oo ^/'^M(/)C;(/')log|T(/)-r(/') 



(210) 



In the above formulas it was implied that C(oo) = 0. One could switch to more 
traditional circular parametrization by mapping the upper half plane inside the unit circle 



1-UJ 

^1 + 1^2 — — — ; ^ — fe^" ; r < 1. 



(211) 



The real axis is mapped at the unit circle. Changing variables in above integral we find 

de „ f 1 1 ■ 



J-TT TT V 1 



rexp(ia — i6) l + exp(— ^6')_ 



Here 



0(^) = / tan- . 



(212) 



(213) 



The last term represents an irrelevant translation of the surface, so it can be dropped. 
The resulting formula for the area reads 



1 l"K rTT 



je'" — e 

or, after integration by parts and inverting parametrization 



(214) 



m = ^ £ dct> £ dct>'C'^{<i>) C'^{ci>') log 



. e{4>) - em 

sm 



(215) 
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Let us now minimize the functional of the boundary parametrization f{t) 

(we shall stick to the upper half plane). The straightforward variation yields 

which duplicates the above diagonality condition ( |204| ). Note that in virtue of this condi- 
tion the normal vector n^{x) is directed towards 92 at the boundary. Explicit formula 
reads 

Let us have a closer look at the remaining nonlinear integral equation (|216|) . In terms 
of inverse parametrization it reads 

r+oo C'(f)C'(f') 

= ^ / df ..f^;,! , (218) 
Introduce the vector set of analytic functions 

F,{z) = / (219) 

which decrease as z^^ at infinity. The discontinuity at the real axis 

QF,{T + tO) = C'^{f)f{T) (220) 

Which provides the implicit equation for the parametrization /(r) 

J dr'^F^{r + tO) = C,{f) (221) 

We see, that the imaginary part points in the tangent direction at the boundary. As for 
the boundary value of the real part of F^{t) it points in the normal direction along the 
surface 

^F^ oc (222) 
Inside the surface there is no direct relation between the derivatives of X^{^) and F^{^). 
The integral equation ( ^16| ) reduces to the trivial boundary condition 

F2(t + zO) = F2(t-^0) (223) 

In other words, there should be no discontinuity of F^ at the real axis. The solution 



compatible with analyticity in the upper half plane and z decrease at infinity is 



F'^iz) = {1 + ur Piu); u = —-^ (224) 



I — z 

I + z 

where P{oo) defined by a series, convergent at |ci;| < 1. In particular this could be a 
polynomial. The coefficients of this series should be found from an algebraic minimization 
problem, which cannot be pursued forward in general case. 
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The flat loops are trivial though. In this case the problem reduces to the conformal 
transformation mapping the loop onto the unit circle. For the unit circle we have simply 

Ci + iC2 = uj; Fi=iF2^ ~^^"T^' ^ " °" ^^^^^ 

Small perturbations around the circle or any other flat loop can be treated in a systematic 
way, by a perturbation theory. 



I Kolmogorov triple correlation and time reversal 

Are there any restrictions on the circulation PDF from the known asymmetry of velocity 
correlations, in particular, the Kolmogorov triple correlation? The answer is that the 
Kolmogorov correlation does not imply the asymmetry of vorticity correlations. 

Taking the tensor version of the | law in arbitrary dimension d 

E ( 2 \ 

(t^a(0)f/3(0)v^(r)) = V "^'^'^ ^ ^^^'^^ ~ ~d^''^^V ^^^^^ 

and differentiating, we find that 

(^a(0)^;M0)c^7A(r-)) = (227) 

So, the odd vorticity correlations could, in fact, be absent, in spite of the asymmetry of 
the velocity distribution. 
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